Analysis of the different interventions scenario for programmatic measles control in Bangladesh: A modelling study

In recent years measles has been one of the most critical public health problem in Bangladesh. Although the Ministry of Health in Bangladesh employs a broad extension of measles control policies, logistical challenges exist, and there is significant doubt regarding the disease burden. Mathematical modelling of measles is considered one of the most effective ways to understand infection transmission and estimate parameters in different countries, such as Bangladesh. In this study, a mathematical modelling framework is presented to explore the dynamics of measles in Bangladesh. We calibrated the model using cumulative measles incidence data from 2000 to 2019. Also, we performed a sensitivity analysis of the model parameters and found that the contact rate had the most significant influence on the basic reproduction number R0. Four hypothetical intervention scenarios were developed and simulated for the period from 2020 to 2035. The results show that the scenario which combines enhanced treatment for exposed and infected population, first and second doses of vaccine is the most effective at rapidly reducing the total number of measles incidence and mortality in Bangladesh. Our findings also suggest that strategies that focus on a single interventions do not dramatically affect the decline in measles incidence cases; instead, those that combine two or more interventions simultaneously are the most effective in decreasing the burden of measles incidence and mortality. In addition, we also evaluated the cost-effectiveness of varying combinations of three basic control strategies including distancing, vaccination and treatment, all within the optimal control framework. Our finding suggested that combines distancing, vaccination and treatment control strategy is the most cost-effective for reducing the burden of measles in Bangladesh. Other strategies can be comprised to measles depending on the availability of funds and policymakers’ choices.


Introduction
Some death-dealing contagious diseases infected many people in different areas of the world in the last few decades. Measles, a human respiratory disease, is one of these devastating The concept of mathematical modeling is used to describe different real-world phenomena's. The researchers have been using mathematical modeling in epidemiology since the middle of the 20th century to interpret the mechanism of spread and find ways of controlling and eliminating infectious diseases. Some renowned researchers have already done a few important studies on measles transmission dynamics deterministically using various essential compartments in recent decades. In Ochoche and Gweryina [16], a SIR deterministic compartmental model of measles was studied, highlighting the impact of vaccination. The authors realized that it is possible to eliminate measles through vaccination, and they advised that at least two doses of vaccine should be mandatory for all children. For controlling the transmission dynamics of measles, an optimal vaccine coverage level was investigated by Momoh et al. [17]. Edward et al. discussed a deterministic mathematical model to eradicate and control the measles pandemic [6]. They observed that the measles would be died out if the basic reproduction number is less than one, i.e. R 0 <1. In another study Fred, Sigey [18], the authors performed on an SEIR mathematical model with vaccination, and they suggested that the mass vaccination and the detection of measles cases as early as possible should be introduced to abet transmission.
Additionally, Christopher et al. studied a mathematical measles model with vaccination and measles drug therapy [19]. They found that the vaccination in susceptible class and drug therapy in exposed class to identified infected person eliminated measles more efficiently. In the study [20], another SEIR mathematical model also investigated to control the transmission dynamics of measles. Although they discussed stability, disease-free and endemic equilibriums, they did not address the impact of the vaccination in their study. In the abovementioned studies, all researchers have developed a deterministic mathematical measles model considering different compartments, including vaccination, to describe the transmission dynamics of measles. Nevertheless, no research has not been done so far focusing on the mathematical model of the mechanism of the spread of measles with double dose vaccination and different intervention scenarios for programmatic measles control.
In this study, we developed a measles model with a double dose of vaccination to describe the transmission dynamics of measles in Bangladesh. The model is calibrated to the Bangladesh demographic and measles incidence data from 2000 to 2019 to estimate the key transmission and progression parameters. Multiple intervention scenarios were considered to explore the impact of each on its own and when combined on measles incidence and mortality. We also performed cost-effective analysis to explore which intervention will be most cost-effective compare to others. This study depicts Bangladesh specific elimination strategies and describes the results of different levels of investment in future on measles control.
The entire paper is designed as follows: Section 2 introduces model development, basic reproduction number estimation, parameter estimation and sensitivity analysis. In Sections 3 and 4, we discussed the development of different scenarios and cost-effective strategy. Finally, section 5 contains discussion and conclusions.

Model development
We developed a compartmental transmission dynamics measles model between the following mutually exclusive compartments: susceptible individuals, S(t); those who have not yet infected with the disease but might become infected; first dose vaccinated individuals, V 1 (t); those who have received the first dose of vaccine; second dose vaccinated individuals, V 2 (t); those who have received the second dose of vaccine; Exposed individuals, E(t); representing those that are infected and have not yet developed active measles disease; Infected individuals, I(t); those who are infected and infectious; and Recovered individuals, R(t); those who were previously infected and successfully recovered.
The total population size N(t) is assumed to be constant and well mixed: To ensure the population size constant, we replace all deaths as newborns in the susceptible compartment. It includes death through natural causes, which occurs in all states at the constant per-capita rate μ, and measles-related deaths, which occur at the constant per-capita rate δ. Susceptible population (S) who receive the first dose of vaccine move to the vaccinated compartment at a rate η. The first dose of vaccinated population V 1 moves to the susceptible compartment at a rate ρ, and the rest of the population moves to the second dose of vaccinated population V 2 at a per-capita rate σ. The second dose of the vaccinated population also moves to the recovery compartment at a rate ω. Individuals in the S compartment may be infected with the measles virus at a rate λ = βSI, where β is the transmission rate between infected and susceptible population. Once infected, individuals move to the exposed compartment E. A proportion of the exposed population progress to the infected compartment at a per-capita rate α, and the rest of the exposed population progress to the recovered compartment at a percapita rate κ due to the treatment of the exposed population. The proportion of the infected individuals move to the recovery compartment due to the treatment rate τ and natural recovery rate γ. The model flow diagram is presented in

PLOS ONE
Scenario analysis for programmatic measles control in Bangladesh: A modelling study From the above-mentioned, the transmission dynamics of measles is given by the following compartmental system of nonlinear ordinary differential equations that describe the model: The initial conditions of the system (2)-(7) are of the form It can be easily shown that the solution of the system (2)-(7) subject to the initial conditions (8) exists and is nonnegative for all t�0.

Basic reproduction number
The basic reproduction number can be determined using the method of next-generation matrix [21]. The next-generation matrix is the production of matrices T and −S −1 where the matrix T represents the rate of infection transmission in E and I compartments and the matrix S describes all other transfer across the compartments. The matrices T and S are given as The next-generation matrix is The basic reproduction number is the Eigen-value of largest magnitude of the next-generation matrix (K). Hence the basic reproduction number is obtained as

Parameter estimation
We estimated the measles model parameters from fitting different combinations of parameters in Eqs (2)- (7) to the actual number of measles cases in Bangladesh from 2000 to 2019 [22]. In order to parameterize measles model (2)-(7), we obtained some of the parameter values from the literature (see Table 1), rest of the parameters were estimated from data fitting (see Fig 2). The estimation of parameters was carried out using the least-squares method, which minimises summation of the square errors given by ∑(A(t,y)−B actual ) 2 subject to the measles model (2)-(7), where B actual is the actual reported measles data, and A(t,y) denotes the solution of the model corresponding to the number of measles cases over time t with the set of estimated parameters, denoted by y.

Sensitivity analysis
We perform the sensitivity of the model basic reproduction number (R 0 ) to the model parameters using the Latin Hypercube Sampling (LHS) method with 10000 runs per simulation. The LHS is a Monte Carlo stratified sampling technique that permits us to concurrently achieve an unbiased assessment of the model output for a particular set of input parameter values. The Partial Rank Correlation Coefficients (PRCCs) for the full range of parameters are shown in Fig 3. The PRCCs for the basic reproduction number in Fig 3 is produced using the expressions R 0 . Results show that parameters transmission rate (β) and progression rate (α) have a positive correlation with the model outcomes R 0 , which means that decreasing these parameters values will reduce the prevalence of measles. On the other hand, parameters ρ, σ, γ, κ, τ and η have a negative correlation with the model outcomes R 0 , which indicates that increasing these parameters will decrease the outbreak of measles.

Ethical approval
This study based on aggregated measles surveillance data in Bangladesh taken from the World Health Organization. No confidential information included because analyses were performed at the aggregate level. Therefore, no ethical approval is required. γ Natural recovery rate 0.6 [2] κ Treatment rate for exposed population 0.08 [26] τ Treatment rate for infected population 0.14 [26]

Scenario development and analysis
This section developed four potential intervention scenarios to explore the dynamics of measles incidence and mortality in Bangladesh. These scenarios are detailed in Table 2, 3, 4 and 5.
We parameterized these proposed responses to our model structure to assess the effect of these responses during the period 2020-2035. In 2014, Bangladesh, one of 11 countries in the South-East Asia Region, adopted a national goal for measles elimination by 2018. In Bangladesh the Ministry of Health achieved its objectives in first and second doses vaccines, effective treatment and overall management through partnership, engaging all care providers (GO-NGOs) and making available free vaccination support. Estimated coverage with the first dose vaccine increased from 74% in 2000 to 94% in 2016. The second dose vaccine was introduced in 2012, and its coverage increased from 35% in 2013 to 93% in 2016 [25]. In Bangladesh, treatment coverage for asymptomatic and symptomatic measles patients are around 8% and 14%, respectively [26]. Therefore, it is essential to increases the treatment rate for both asymptomatic and symptomatic measles patients in Bangladesh. Scenario 1 simulates a continuation of the programmatic situation during the period 2020-2035. During this time, we incorporated four different intervention strategies: increasing both first and second dose vaccine; improving treatment rates for exposed (asymptomatic) and infected (symptomatic) population. We implement these as single interventions and compare

PLOS ONE
Scenario analysis for programmatic measles control in Bangladesh: A modelling study them with baseline (see Table 2) to explore the impact of each intervention on total measles incidence and mortality. In this scenario, the first and second dose vaccine increase from 94% and 93% to 100% and 100%, respectively. Further, treatment rates for asymptomatic and symptomatic cases improved through a combination of case-finding strategies and improved knowledge of standard operating procedures for measles treatment commencement. The case finding is to identify of symptomatic patients attending a health facility, either their initiative or referred by another health facility, health worker, and community volunteer. These strategies improve the treatment rate for an asymptomatic case from around 8% to 100% and symptomatic cases from around 14% to 100%. Scenario 2 consists of a combination of six-double intervention strategies, which include (i) enhanced first dose and second dose vaccine, (ii) enhanced the first dose vaccine and treatment rate for exposed population, (iii) enhanced first dose vaccine and treatment rate for infected population, (iv) improved second dose vaccine and treatment rate for an exposed population, (v) improved second dose vaccine and treatment rate for an infected population, and (vi) enhanced treatment rate for exposed and infected population.
Scenario 3 incorporates a combination of four-triple intervention strategies such as (i) first and second dose vaccine as well as treatment rate for an exposed population, (ii) first and second dose vaccine as well as treatment rate for an infected population, (iii) first dose vaccine, treatment rate for an exposed and infected population, (iv) second dose vaccine, treatment https://doi.org/10.1371/journal.pone.0283082.g003 rate for exposed and infected population. Finally, scenario 4 describes the broad scale-up combination of first and second dose vaccines and the treatment rate for exposed and infected populations. Under this scenario, all model parameters values changed simultaneously and considered overall scale-up values for each parameter. However, each category of intervention could involve several potential specific activities. For example, treatment for exposed and infected populations could include training of doctors, nurses, and pharmacists on measles guidelines, monitoring and managing supplies of high-quality drugs. We assumed that these different activities would have a similar impact on the model and did not model the impact of these explicit activities independently.
The estimated outcomes in this study, over a 15-year time frame, included: first dose and second dose vaccine, treatment for exposed and infected population. To better understand the contribution of changing specific interventions with different scenarios, estimated outcomes were presented separately for each intervention assumed to be influenced by the scenario. Results from scenario one are presented in Figs 4 and 5 as well as Table 2. From scenario 1, we observed that within a four-single intervention strategy, the treatment rate for the exposed population is the most effective than other single interventions, which reduces more measles incidence and mortality (see Figs 4(A3) and 5(B3), and Table 2) in Bangladesh. Hence, it is the preferable single intervention strategy. Alternatively, the treatment rate for the infected population is another option. Figs 6 and 7, Table 3, represent scenario 2, which includes our proposed double interventions strategies. Each of the interventions resulted in decreasing the number of measles incidence and mortality. The analysis shows that a combination treatment rates of the exposed and infected population are the best dual intervention strategy for reducing the number of measles incidence and mortality in Bangladesh (see Figs 6(A6) and 7(B6), and Table 3). Alternative, the combination of first dose vaccination rate and treatment rate for an exposed population is another option. Figs 8 and 9 and Table 4 depict scenario 3, which include a combination of triple interventions strategies. As expected, each of the strategies decreased the number of measles incidence and mortality in Bangladesh. The combination of first dose vaccination rate, treatment rates for the exposed and infected population is the best triple intervention strategy (see Figs 8(A3) and 9(B3), and Table 4). Alternative, second dose vaccination rate, and treatment rates for exposed and infected populations are also better triple intervention strategies.

PLOS ONE
Scenario analysis for programmatic measles control in Bangladesh: A modelling study Finally, Fig 10 and Table 5 represents scenario 4, which includes first and second dose vaccination rates and treatment rates for exposed and infected populations. Under scenario 4, measles incidence and mortality reduce enormously over the 15-years due to the combination of immensely expanding first and second dose vaccination rates and treatment rates for exposed and infected populations. We also have compared all the scenarios to know which is the most effective (see Table 5). Our finding suggests that within the best scenarios analysis, scenario 4 is the most effective, which reduce the massive number of measles incidence and mortality in Bangladesh. However, depending on the availability of fund, other scenario in Table 5 can be considered.

Optimal control strategy and cost-effective analysis
In this section, we performed optimal control strategy and implemented three time dependent control variables to explore their effectiveness and cost-effective analysis in controlling the spread of measles in the population of Bangladesh. The time dependent control variables u 1 (t), u 2 (2) and u 3 (t) are defined as follows: i. u 1 (t) denotes the distancing control strategy that is the effort at inhibiting the virus transmission from Exposed and infected population. This can be reached through public health advocacy for social distancing, good personal hygiene, diagnosis campaigns, and education programs for public health. Noting that u 1 (t) = 1 indicates the policy effectively protects against infection, while u 1 (t) = 0 denotes strategy failure.

PLOS ONE
Scenario analysis for programmatic measles control in Bangladesh: A modelling study ii. u 2 (t) represents the vaccination control strategy, it is assumed that the number of vaccines available during this time period and they are all administrated and used completely. If u 2 (t) = 1, then the control strategy is effectively used, while u 2 (t) = 0 means the absence of the control strategy.
iii. u 3 (t) indicates control variable to enhance the treatment of infected population with a view to ensure the rapid provision of additional treatment includes providing comfort measures to relieve symptoms and preventing complications. Perceiving u 3 (t) = 1, then the control strategy is effectively treating the disease, while u 3 (t) = 0 means the strategy failure.
Subsequently, the optimal control model with the three above-mentioned time-dependent variables is given by the following non-linear differential equations: The goal of presenting the three control variables is to seek the optimal solution required to minimise the numbers of Exposure and infected individuals at minimum cost. Hence, the objective function for this optimal control problem is given by where, constants A i , i = 1,2 are positive weights essential to balance the objective function. Following other works on infectious diseases control problem [27][28][29][30], quadratic cost on the controls are chosen to ensure the control has only one extremum (i.e. maximum or minimum), where 1 2 B 1 u 2 1 t ð Þ is the total cost of executing the distancing, and 1 2 B 2 u 2 2 t ð Þ is the total cost of vaccination and 1 2 B 3 u 2 3 t ð Þ is the total cost of treatment for infected individuals over the time interval [T 0 , T f ] (where initial time T 0 = 0, final time T f = 15 years period).  and σ, (A2) η and κ, (A3) η and τ, (A4) σ  and κ, (A5) σ and τ, and (A6) κ and τ. https://doi.org/10.1371/journal.pone.0283082.g006 Precisely, the optimal control strategy u * ¼ ðu where, U is the non-empty control set defined by U ¼ fðu 1 ; u 2 ; u 3 Þ : ðu 1 ðtÞ; u 2 ðtÞ; u 3 ðtÞÞ are measurable with 0 � u 1 ; u 2 ; u 3 � 1 for t 2 ½T 0 ; T f �g: .

PLOS ONE
Scenario analysis for programmatic measles control in Bangladesh: A modelling study the model (9). The expected outcome for minimising control problem as performed in [30,32] is adapted below. Now using Pontryagin's maximum principle, we obtain the following theorem. Theorem: Given that ðu * 1 ; u * 2 ; u * 3 Þ minimises the objective function (10) subject to the corresponding system (9), then the adjoint variables λ i , i = 1,2,3,. . .,6, satisfy the following system with the terminal (transversality) conditions Further, the optimal control pair ðu * 1 ; u * 2 ; u * 3 Þ is given as follows Proof: The existence of the optimal controls u * 1 ; u * 2 and u * 3 such that Jðu * 1 ðtÞ; u * 2 ðtÞ; u * 3 ðtÞÞ ¼ min U Jðu 1 ; u 2 ; u 3 Þ with state system (9) is given by the convexity of the objective function integrand. By Pontryagin's Maximum Principle [31], the adjoint equations and transversality conditions are obtained. Differentiation of Hamiltonian H 1 for the state variables gives the following system, with λ i = 0, for i = 1,2,3,. . .,6.
Optimal controls u * 1 ðtÞ; u * 2 ðtÞ and u * 3 ðtÞ are derived by the following optimality conditions, at u * 1 ðtÞ; u * 2 ðtÞ and u * 3 ðtÞ on the set U. On this set This ends of the proof.
Here, we implemented Runge-Kutta fourth-order forward and backward method using MATLAB programming language to solve the subsequent optimality system which consists of (9) and (13) with the characterization (15) within the period of [0, 15] years. The weight constants adopted for balancing the objective function (10) are selected to ensure that no term dominates the other. Therefore, we used equal weight constant for minimising the infectious classes, so that A 1 = A 2 = 1. Under other conditions, the weight constants for determining efforts or cost essential to implement the controls are comparatively different, and outcomes in values for B 1 = 50, B 2 = 100 and B 3 = 150 are consistent with previous modelling research. [33]. Details of the numerical procedure for simulating the obtained optimality system are contained [34].
Figs 11 and 12 establish how distancing control, u 1 (t) and vaccination control strategies, u 2 (t) affect the spread of the measles in Bangladesh. As shown in Figs 11 and 12, to minimise the objective function (10), the optimal control u 1 (t) and u 2 (t) are continued at the maximum level (i.e. 100%) for about 8 years and 5 years respectively for Bangladesh population before relaxing to the minimum in the final time. Also as expected, the number of measles infectious individuals are reduced when control is in place. On the other hand, Fig 13 displays the effects of treatment control strategy u 3 (t) on the dynamics of measles infection in Bangladesh. We observed that treatment control strategy has small impact on exposed and infected population while distancing control strategy has high impact for reducing the burden of measles cases. Fig 14 shows the implication of combining the three optimal controls in bringing down the total number of infectious human to zero in Bangladesh. It is observed that optimal solution has achieved when distancing control strategy (u 1 ) is strictly followed to at the maximum level of 100% for 7 years, while the vaccination and treatment control strategies (u 2 , u 3 ) are at a maximum level above 60% and 20% respectively. It can be seen that the combination of the three control strategies is significantly more effective to decrease the spread of the measles compare to implement each control strategy individually, which is consistent with the previous modelling studies [30,35,36].

PLOS ONE
It is crucial to identify the most cost-effective strategy among distancing, vaccination and treatment control as well as their combination control strategies to optimally mitigate the spread of measles at the possible minimum cost. This is performed by associating the differences among each intervention's costs and outcomes; obtained by estimating the incremental costeffective ratio (ICER), which is defined as the extra cost per additional intervention outcome. Incrementally, when analysing two or more competing intervention policies, one intervention is associated with the next less effective option. The ICER numerator is given by the total difference in intervention costs, active measles cases averted costs and averted productivity losses if applicable, between each scenario and baseline. The ICER denominator is the difference in total number of active measles cases averted. Hence, the following formula obtains the ICER:

ICER ¼
Difference in total costs between control strategies Difference in total number of active cases averted by control strategies ; ð16Þ

PLOS ONE
We also performed the average cost-effectiveness ratio (ACER), which estimates the effectiveness of a particular intervention technique. The ACER is the ratio between the total cost incurred and the total number of active measles cases averted by that strategy. This is calculated by The total cost for each of distancing, vaccination and treatment implementation and mutual effort of the optimal control strategy is obtainable from the objective function (10). The cases averted is invaded by computing the difference between infectious individuals with and without control strategy. Let, P 1 , P 2 , P 3 and P 123 respectively represent single distancing control strategy u 1 (t), single vaccination control strategy u 2 (t), single treatment control strategy u 3 (t) and the combined effort of the three strategies. Table 6 summarises the ICER and ACER for each and the combination of the control variables u 1 (t), u 2 (t) and u 3 (t) in increasing order of the total infection averted. The ICER and ACER results for P 1 , P 2 , P 3 and P 123 are calculated using (16) and shown in Table 6 follows.
Comparing P 1 , P 2 , P 3 and P 123 in Table 6 and Fig 15, it is seen that combined control strategy P 123 is most cost-effective which reduce a significant number of measles cases with low cost compared to P 1 , P 2 and P 3 individually, while P 3 is the least cost-effective intervention strategy among them.

Discussion and conclusion
Bangladesh is a resource-poor, high burden measles country compared to other South-East Asian countries, and the transmission dynamics and epidemiology of measles are poorly understood. Therefore, due to paying less attention to the spread of measles, Bangladesh has been suffering from measles for many years even though measles is a vaccine-preventable disease with adequate and timely vaccination [37]. Indeed, ensuring vaccines for all people is a high burden and expensive for Bangladesh, as in other South Asian countries, but still the collected vaccines are not administered properly and therefore a considerable number of vaccines are wasted [38]. Even at the beginning of the last decade, more than half (about 69.7%) of measles vaccines distributed at both ward and village levels almost in all districts were wasted due to improper distribution planning or supply process and lack of proper maintenance [39].
Owing to these problems, the measles vaccination rate was comparatively low than that of other diseases vaccinations [40]. However, subsequently, the Bangladesh government set a target of achieving more than 95% measles vaccination coverage by 2018 to eliminate measles [41]. Finally, following WHO measles elimination strategies, Bangladesh has achieved around 93% coverage with second vaccine doses in every district and delivered through the routine immunization program and financial support from the treatment program and programmatic strength [25]. Although measles control in Bangladesh has significantly progressed-improved first and second doses of vaccine, availability of free treatment services, involvement of multiple partners, newer diagnostic facilities, sufficient human resources, adequate capacity and guidelines-more efforts are required.
In this paper, we presented a measles model with first and second dose vaccine to explore the transmission dynamics of measles in Bangladesh. We derived the basic reproduction number of the measles model and found that basic reproduction numbers play an essential role in the dynamics of measles outbreaks. We fitted measles incidence data with the model from the WHO report to estimate the model parameters. Sensitivity analyses were performed to determine the relative importance of several parameters used in our model. It will help epidemiologists and public health officials focus on the more significant parameters for formulating the measles control strategy. Our analysis led to the observations that transmission rate significantly contributes to the basic reproduction numbers of measles dynamics and needs to be calculated precisely for an accurate outcome.
We acknowledged the importance of comprehensive countrywide programmatic improvements to measles control in Bangladesh. Without such extensive approaches, further rises in the overall disease burden are expected. Four different scenarios were considered (scenario 1, scenario 2, scenario 3 and scenario 4) from combinations of first and second dose vaccines and treatment for exposed and infected populations. Numerous settings were examined to measure the effectiveness of the scenario strategies. In scenario 1, within the four single-intervention strategies, the treatment for the exposed population is the most effective intervention compared to other single-intervention strategies. The result shows that if extensive programmatic development does not occur, especially with the expansion of treatment rate for the exposed population, it is likely that measles will remain uncontrolled in Bangladesh. However, the first and second dose vaccination rates already have covered about 94% and 93%, respectively, in Bangladesh. Therefore, increasing the treatment rate for exposed populations is the better option for controlling measles outbreaks.
Within the six dual-intervention strategies (scenario 2) with treatment for exposed and infected populations is the most effective intervention strategy that reduces the number of incidence and mortality in Bangladesh. Alternatively, activities that aim to improve the first dose of vaccine and exposed or infected population treatment rate are other effective options. However, when we scaled up the treatment rate for the exposed and infected population under scenario 2, the reduction in measles incidence and mortality achieved. Therefore, our outcomes suggest that an improved treatment rate for the exposed and infected population would be insufficient to reduce measles incidence and mortality.
Considering the four-triple intervention strategies (scenario 3), improving the treatment rate for exposed and infected populations with the first dose of vaccination rate is the most effective at measles incidence and mortality reduction compared to other triple intervention strategies. From the analysis of all the scenarios, we found that scenario 4 (first and second dose vaccination rate and treatment rate for exposed and infected population) is the most effective intervention strategy and followed by scenario 3 (first dose vaccination rate and treatment rate for exposed and infected population), scenario 2 (treatment rate for exposed and infected population), and scenario 1 (treatment rate for exposed population). Our findings suggest that the treatment rate for the exposed population is the most critical intervention for decreasing measles incidence and mortality in Bangladesh. However, focus on the first dose and second dose of vaccination rate or treatment rate for the infected population alone will not dramatically affect the decline in measles incidence and mortality in Bangladesh. Taking two or more key interventions simultaneously is the most effective way to reduce measles incidence and mortality in Bangladesh.
Scenarios analysis has been applied in different low and middle income endemic countries to control infectious diseases (e.g. measles, tuberculosis, COVID-19) epidemic [42][43][44][45]. Previous studies show that focus on a single intervention strategy will not dramatically affect the decline in infectious disease outbreak but combined two or more interventions is the most effective for reducing the burden of infectious diseases [43,46], which is found to be consistent with our results.Our scenarios describe a variety of potential responses, extending from inaction to extremely ambitious multifactorial strategies. Despite the challenges of delivering effective programmatic measles control in Bangladesh, we believe it is essential to consider such responses. It is because significant impacts from simple public health interventions have previously been demonstrated in resource-limited settings such as Bangladesh. However, the World Health Organization does not presently recommend comprehensive approaches. Our modelling suggests that the high burden of measles incidence and mortality in Bangladesh is likely to increase with the existing programmatic intervention strategy.
We also implemented an optimal control approach via Pontryagin's Maximal principle [31] and formulated the optimal strategies for controlling the measles epidemic in Bangladesh. Three different control strategies were considered including distancing control strategy (u 1 ), vaccination control strategy (u 2 ) and treatment control strategy (u 3 ). Different settings were examined to measure the cost-effectiveness of the control strategies. Between the three-single control strategies, the distancing control strategy (u 1 ) is better in cost-effectiveness than the vaccination and treatment control strategies which reduce a significant number of measles cases in Bangladesh. Therefore, our results suggest that the Bangladesh government should improve distancing control interventions when only one control strategy is used. Naturally, this strategy actively decreases and/or stops the contact between susceptible and infectious individuals of measles. However, combined implementation of distancing, vaccination and treatment strategies is the most cost-effective measure for reducing the burden of measles in Bangladesh.
Optimal control strategies have been applied in other endemic settings to minimise infectious diseases (e.g. measles, tuberculosis, COVID-19) cases and intervention implementation costs. Previous studies show that distancing strategy is the best strategy for the single control strategy implementation to decrease disease burden and intervention costs [33,42,47], which is similar to our results. However, our finding also suggests that combining three control strategy is the most effective way to decrease the measles burden of Bangladesh, consistent with previous works [33,48,49].
In Bangladesh, infectious disease surveillance is not fully recognised and the risk of bias cannot be prohibited. More precise data should be put in place to address alarms related to measles. Precise data leads to better estimation of vital parameters, and this means our projected intervention to decision support is data-dependent. Hence, local and national level policy-makers need to adjust the possibility of under-reporting bias when investigating our results.